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Abstract.  The  shot  noise  in  correlated  single-electron  tunneling  can  be  calculated  using  the 
Langevin  approach.  It  is  equivalent  to  the  Fokker-Plank-type  approach  in  the  “orthodox”  frame- 
work, however,  the  advantage  is  a natural  possibility  to  describe  also  the  fluctuations  in  the  quantum 
frequency  range. 


Correlated  single-electron  tunneling  [ ] remains  an  attractive  topic  during  last  decade. 
The  theoretical  study  of  the  shot  noise  in  this  regime  is  interesting  because  the  noise  reflects 
rather  directly  the  correlations  in  the  charge  transport,  and  also  because  the  noise  determines 
the  ultimate  sensitivity  of  single-electron  devices. 

The  shot  noise  in  single-electron  tunneling  is  due  to  the  randomness  of  tunneling  events. 
The  earlier  developed  formalism  [ | of  noise  calculation  is  of  the  Fokker-Plank  type  and  is 
based  on  the  deterministic  master  equation  of  “orthodox”  theory  [ ] (we  assume  sufficiently 
large  tunnel  resistances,  Rj  ^ h/e^): 


& = Fa , 
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where  the  element  of  the  vector  a is  the  probability  to  find  the  system  in  the  charge  state 
n = {n\,  ...iil]  (which  is  characterized  by  the  numbers  n,-  of  excess  electrons  in  each  of  L 
internal  nodes  of  the  system)  and 
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where  T^^n  are  the  corresponding  tunneling  rates  and  the  summation  over  the  junction 
number  / is  necessary  when  an  electron  can  tunnel  to  (from)  an  internal  node  from  (to) 
different  external  electrodes. 

To  find  the  mutual  spectral  density  for  two  processes  X (t)  and  Y(t)we  can  calculate  first 
the  correlation  function  Kxy  (t)  = {X{t+r)Y  (f))  — (X)  (T)  (brackets  denote  the  averaging 
over  time)  and  then  take  the  Fourier  transform  Sxy((o)  — 2 KxY(t)  exp(iia;r)  dr.  If 
both  X and  Y are  functions  of  the  charge  state  n (for  example,  potential  of  a node)  then  the 
correlation  function  is  given  by  the  simple  expression 

Kxvir)  — Ojr)  X{m)  a{r,m\n)Y(n)cr^‘  +0(-r)  E Y (m)  cr(m,  ^r\n) 

m,n  m,n 

xX(uK'-(X)(T),  (3) 

where  cr(r,  m\n)  is  the  retarded  Green’s  function  of  Eq.  (1)  being  the  probability  to  find 
the  system  in  the  state  m at  f = r > 0 if  at  I = 0 it  was  in  the  state  n,  {X)  — ^ , 

and  is  the  stationary  distribution,  T — 0,  ^ • 
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However,  if  X and/or  Y represent  the  current  through  a tunnel  junction  or  in  an  ex- 
ternal lead,  Eq.  (3)  should  be  modified.  For  example,  if  X (t)  is  the  current  contribution 
corresponding  to  tunneling  events  while  Y {t)  corresponds  to  then  (similar 

to  Ref.  [ ]) 


/^XF(r)/c±c±  = e{x)  a(x,n\m')  +d(-x) 

xcrin, -x\m)ri^n  <^n  - "t"  ^m}n'^nn 


m <—n 


.(4) 


Here  the  last  term  is  responsible  for  the  high-frequency  limit.  The  effective  charges  and 

are  determined  by  the  direction  of  electron  tunneling,  = —e^_,  and  by  the  circuit 
capacitances  [ ] (so  that  — e only  if  the  current  through  junction  / is  measured).  Any 
current-current  correlation  function  can  be  written  as  a sum  of  A'xy(r)  given  by  Eq.  (4) 
over  all  possible  transitions  between  charge  states  (such  a sum  is  a counterpart  of  Eq.  (3) 
in  which  the  sum  is  written  explicitly). 

The  expressions  for  spectral  densities  directly  follow  from  Eqs.  (3)  and  (4)  because 
the  Fourier  transformation  affects  only  the  evolution  operator  a{x,  m\n),  and  the  corre- 
sponding Green’s  function  in  the  frequency  representation  is  simply  obtained  from  Eq.  (1): 
<y{a>,  m\n)  — — E)“'  \mn,  where  1 is  the  unity  matrix.  For  example,  Eq.  (4)  leads 

to  the  following  spectral  density: 
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This  method  allows  to  calculate  all  spectral  densities  within  the  framework  of  “orthodox” 
theory,  and  at  least  for  the  single-electron  transistor  the  numerical  procedure  is  rather  trivial 
[ ] because  the  matrix  E is  three-diagonal  and  the  matrix  inversion  is  straitforward. 

The  alternative  method  of  calculations  can  be  based  on  the  Langevin-type  approach 
[ ].  Because  of  the  discrete  nature  of  the  charge  states,  the  random  term  cannot  be  simply 
added  (in  analogy  with  the  standard  Langevin  method)  to  some  evolution  equation  for  the 
“coordinate”  n(t),  but  should  be  introduced  into  the  master  equation  (similar  to  Ref.  [ ]). 
The  derivation  of  the  formalism  can  be  understood  in  the  following  way.  Let  us  consider 
the  ensemble  of  M (M  ^ 1 ) independent  similar  circuits,  and  let  us  average  all  magnitudes 
over  this  ensemble.  Then  the  average  (over  time)  currents  and  voltages  will  not  change 
(due  to  ergodicity),  but  the  spectral  densities  of  fluctuations  (second  order  magnitudes)  will 
decrease  M times.  Hence,  to  calculate  spectral  densities  of  the  initial  system,  we  can  take 
the  leading  (~  M“’ , M oo)  order  of  the  spectral  density  of  magnitudes  averaged  over 
the  ensemble. 

In  contrast  to  the  single  system,  the  dynamics  of  the  large  ensemble  is  easily  described 
using  the  Langevin  approach.  At  any  moment  of  time  t the  ensemble  can  be  charac- 
terized by  “coordinates”  Ma„{t)  which  represent  the  numbers  of  participants  being  in 
different  charge  states  n (notice  that  now  a„  (t)  is  not  a probability  but  the  fluctuating  co- 
ordinate). While  in  the  stationary  state  the  average  number  of  transitions  from  state  m 
to  state  n during  small  time  At  is  given  by  MFm^nO'n  the  r.m.s.  of  this  number  is 
obviously  (Mr,„^„cr^‘  At)^^^.  Hence,  the  recipe  is  the  following  [ ]:  for  each  average 
flux  in  the  space  of  charge  states,  we  should  add  in  the  master  equation  the 
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random  5-correlated  (white)  flux  with  the  corresponding  “seed”  spectral  density 

given  by  Schottky-like  formula, 

bm  (r ) = ^ ^ "b  (j)  — ^ ^ Hm  (t)  j (6) 

n n,j 


S J J,  (w)  = 2M- 


* m-<—n 


(7) 


For  the  fluxes  in  opposite  directions  (m  ■(—  n and  n m)  we  should  apply  f (r)  for  each 
direction,  so  that  the  random  flux  does  not  vanish  even  if  the  net  average  flux  is  zero. 

Because  of  the  linearity  of  Eqs.  (6)-(7)  the  final  spectral  densities  of  the  averaged  (over 
M)  magnitudes  are  obviously  proportional  to  1 /M.  Hence,  rescaling  to  the  single  system 
can  be  done  formally  assuming  M — \ \n  Eqs.  (6)-(7).  So,  instead  of  keeping  M and 
rescaling  at  the  final  stage,  we  will  use  M — \ \n  all  equation  below. 

Using  the  standard  procedure  we  find  the  Fourier  transform 


Cm  (<W)  = 


(— /o)l 


r)-'l  §„(a>). 

Jm/i 


(8) 


Then  for  the  occupation-occupation  spectral  density  we  obtain  the  expression 


-^cr,ncr/i 


= ^ _ r)-'l  , \iicoi  - r)-'l  , 

L Jinm'  L Jnn' 
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which  coincides  with  the  result  of  Fokker-Plank  approach  (Fourier  transform  of  Eq.  (3) 
without  X and  Y factors). 

The  technique  is  similar  for  the  current-current  fluctuations.  The  case  of  Eqs.  (4)  and 
(5)  corresponds  to  currents 


X{t) 


,an it)  + (o]  , Y(t)  = ~e{  [r£,^„,cr„, it)  + it) 


, (10) 


and  the  straitforward  calculations  using  Eqs.  (7)  and  (8)  lead  to  Eq.  (5).  Thus,  the  Fokker- 
Plank  method  is  equivalent  to  the  Langevin  method  within  the  “orthodox”  framework. 
However,  in  contrast  to  the  former  approach,  the  Langevin  method  easily  allows  phe- 
nomenological generalization  for  the  fluctuations  in  the  range  of  high  (“quantum”)  fre- 
quencies, hw  ^ eV . 

Let  us  remind  that  in  “orthodox”  theory  [ ] the  tunneling  rate  F = /o(tF/e)/e[l  ^ 
exp(^tT/7’)|  is  determined  by  the  energy  gain  due  to  tunneling  W — eVh  — jlCeff, 
where  loiv)  is  the  “seed”  I-V  curve  of  the  junction  (in  the  linear  case  loiv)  — v/Rj), 
Vb  is  the  voltage  across  the  junction  before  the  tunneling,  Ceff  is  the  effective  junction 
capacitance  (which  also  accounts  for  the  environment),  and  T is  the  temperature.  The 
generalization  of  the  Langevin  method  is  the  substitution  of  Eq.  (7)  by  the  equation 

S j i'  i(a)  — 8mm'8„„'Sjji[r~^  + r~]cr^‘  , 

iQji'^i^n/e  ± Yuo/e) 
e 1 - exp  (-iwi^n  ± hco)/ T ^ 


(11) 
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which  is  derived  for  the  individual  tunneling  event  within  the  standard  tunneling  hamiltonian 
technique  averaging  the  quantum  current-current  correlator  and  then  taking  the  Fourier 
transform.  Eq.  (11)  can  be  considered  as  a generalization  of  the  fluctuation-dissipation 
theorem  and  equations  of  Ref.  [ ] for  the  case  of  single-electron  tunneling.  (Actually,  the 
only  difference  is  that  we  separate  fluctuations  corresponding  to  two  directions  of  tunneling. 
In  absence  of  the  Coulomb  blockade  they  could  be  summed  together  leading  to  the  standard 
factor  [ 1 coth((IF  ± hco)/2T)  instead  of  the  denominator  of  Eq.  (1 1).) 

Equations  (6)  and  (11)  represent  a phenomenological  generalization  in  which  the  low 
frequency  behavior  is  treated  by  the  master  equation  while  high  frequency  properties  are 
taken  into  account  for  individual  tunneling  events.  At  high  frequencies,  o)  ^ E,  the 
occupation-occupation  and  occupation-current  spectral  densities  vanish,  while  for  the 

current-current  spectral  density  instead  of  Eq.  (5)  we  get  j ,/  (a>) , 

because  the  first  terms  of  Eq.  ( 1 0)  are  too  slow  to  give  a contribution.  This  result  coincides 
with  the  result  of  Ref.  [ |.  The  advantage  of  the  Langevin  approach  is  the  possibility  to 
obtain  spectral  densities  in  the  “orthodox”  and  “quantum”  frequency  ranges  using  the  same 
formalism  while  before  they  were  necessarily  being  treated  on  different  footing. 

Notice  that  in  the  “quantum”  frequency  range  the  current  spectral  density  does  not 
correspond  directly  to  the  available  power  because  of  the  contribution  from  zero-point 
oscillations.  The  contribution  to  be  subtracted  from  Su  is  equal  to  2fi,a)ReG(o)),  where 

sm*-n 

the  active  conductance 

ReG(<w)  = ((4)V2fiw)  |f  + - f-|  (12) 

corresponds  to  the  lowest  order  of  photon-assisted  tunneling. 
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